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1. Introduction 

It is already a common knowledge that the accuracy of data for DIS structure functions (SFs) 
allows one to study <2 2 -dependence of logarithmic QCD-inspired corrections and those of power- 
like (non-perturbative) nature independently (see for instance [^] and references therein). And this 
aspect is crucial for the analysis to be performed within some well defined scheme. 

In this contribution we present the results of our recent analysis [JTJ| of DIS SF F 2 (x, Q 2 ) carried 
out over SLAC, NMC, BCDMS and BFP experimental data ^ at NNLO of massless perturbative 
QCD. As in our previous papers [Q, ||| ||] the function F 2 (x, Q 2 ) is represented as a sum of the 
leading twist F 2 ® CD (x, Q 2 ) and the twist four terms: 

F 2 (x,Q 2 )=Ff CD ( X ,Q 2 )(l + ~^-). (1.1) 



Q 1 

As is known there are at least two ways to perform QCD analysis over DIS data: the first one 
(see e.g. ^ ||]) deals with Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) integro-differential 
equations and let the data be examined directly, whereas the second one involves the SF mo- 
ments and allows performing an analysis in analytic form as opposed to the former option. In this 
work we take on the way in-between these two latter, i.e. analysis is carried out over the moments 
of SF F* (x, Q 2 ) defined as follows 

M p QC D/twm/... , q2) = r'^pPQco/twisty... , Ql)dx (L2) 
Jo 

and then reconstruct SF for each Q 2 by using the Jacobi polynomial expansion method 10]. The 
theoretical input can be found in the papers []6, 11]. 



2. A fitting procedure 

The fitting procedure largely follows that used in With the QCD expressions for the Mellin 
moments M^(Q 2 ) analytically calculated according to the formula? given above the SF F 2 (x,Q 2 ) 
is reconstructed by using the Jacobi polynomial expansion method: 

Nmax n 

F 2 k (x,Q 2 )=x«(l -xf £ 0f (*) £cf\a,l5)M* +2 (Q 2 ), 

n=Q j=Q 

where &"' b are the Jacobi polynomials and a, b are their parameters to be fitted. A condition im- 
posed on the latter is the requirement of the error minimization while reconstructing the structure 
functions. 

Since a twist expansion starts to be applicable only above Q 2 ~ 1 GeV 2 the cut Q 2 > 1 GeV 2 
is imposed on the experimental data throughout. The MINUIT program [|12|] is used to minimize 
two variables 

2 o ex p - iy eor 

Xslope — AQexp 

where D = d In F2/0! In In Q 2 . The quality of the fits is characterized by # 2 /DOF for the SF F 2 . 
Analysis is also performed for the slope D that serves the purpose of checking the properties of fits. 

We use free normalizations of the data for different experiments. For a reference set, the most 
stable deuterium BCDMS data at the value of the beam initial energy E = 200 GeV is used. 
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Q 2 


Nof 
points 


HTC 


* 2 (F 2 )/DOF 


a, (90 GeV 2 ) ± stat 


«v(Ml) 


1.0 


797 


No 


2.20 


0.1767 ±0.0008 


0.1164 


2.0 


772 


No 


1.14 


0.1760 ±0.0007 


0.1162 


3.0 


745 


No 


0.97 


0.1788 ±0.0008 


0.1173 


4.0 


723 


No 


0.92 


0.1789 ±0.0009 


0.1174 


5.0 


703 


No 


0.92 


0.1793 ±0.0010 


0.1176 


6.0 


677 


No 


0.92 


0.1793 ±0.0012 


0.1176 


7.0 


650 


No 


0.92 


0.1782 ±0.0015 


0.1171 


8.0 


632 


No 


0.93 


0.1773 ±0.0018 


0.1167 


9.0 


613 


No 


0.93 


0.1764 ±0.0022 


0.1163 


1.0 


797 


Yes 


0.98 


0.1772 ±0.0027 


0.1167 



Table 1: a s (M|) and % 2 in the case of the combined analysis (HTCs stands for higher twist corrections). 

3. Results 

Since the gluon distribution function is not taken into account in the nonsinglet approximation, 
the analysis is substantially easier to conduct; hence the cut on the Bjorken variable (x > 0.25) 
imposed where gluon density is believed to be negligible. The starting point of the evolution is 
taken to be <2o = 90 GeV 2 . These <2q values are close to the average values of Q 2 spanning the 
corresponding data. The previous experience tells us that the maximal value of the number of 
moments to be accounted for is N majc = 8 [Q] (though we check N max dependence just like in the 
NLO analysis) and the cut 0.25 < x < 0.8 is imposed everywhere. 

In [|6|, [[]] the cuts on the kinematic variable Y = (Eq —E)/Eq have been imposed so as to 
exclude BCDMS data with large systematic errors. Here Eq and E are lepton initial and final 
energies, respectively. Upon excluding the set of data with large systematic errors considerably 
higher values of a s (M|) are obtained and rather mild dependence of its values on the choice of Y 
cut is observed. For more details we refer to [§] g|. Once these cuts are applied, a full set of data 
consists of 797 points. 

To verify a range of applicability of perturbative QCD we start with analyzing the data without 
a contribution of twist-four terms (which means F2 = F 2 P ® CD ) and perform several fits with the cut 
Q > Qmin gradually increased. From Table 1 it is seen that unlike the NLO analysis the quality 
of the fits starts to appear fairly good from Q 2 = 3 GeV 2 onwards (at NLO, it starts at Q 2 = 10 
GeV 2 ^) . Then, the twist-four corrections are added and the data with the usual cut Q 2 > 1 GeV 2 
imposed upon is fitted. It is clearly seen that as in the NLO case (see [|]]) here the higher twists do 
sizably improve the quality of the fit, with insignificant discrepancy in the values of the coupling 
constant to be quoted below. 

The parameter values of the twist-four term are presented in Table 2. Note that these for H2 
and D2 targets are obtained in separate fits by analyzing SLAC, NMC and BCDMS datasets taken 
together. For illustrative purposes we visualize those for the hydrogen data in Fig. 1, where the 
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X 


hi, (x) for H2 ± stat 


h^{x) for D2 ± stat 


0.275 


-0.183 ±0.020 


-0.197 ±0.009 


0.35 


-0.149 ±0.028 


-0.171 ±0.015 


0.45 


-0.182 ±0.029 


-0.033 ±0.031 


0.55 


-0.236 ± 0.052 


0.142 ±0.057 


0.65 


-0.180 ±0.135 


0.295 ±0.108 


0.75 


-0.177 ±0.182 


0.303 ±0.158 



Table 2: HTC parameter values obtained in NNLO analysis. 
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Figure 1: Comparison of the HTC parameter h^(x) obtained at LO, NLO and NNLO for hydrogen data (the 
bars indicate statistical errors). 

HTCs obtained at NLO and NNLO levels are seen to be compatible with each other within errors. 

We would like to note that the cut of the BCDMS data, which has increased the a s values (see 
Fig. 1 in [|l|]) improves considerably agreement between perturbative QCD and experimental data. 
Indeed, the HTCs, that are nothing else but the difference between the twist-two approximation (i.e. 
pure perturbative QCD contribution) and the experimental data, are seen to become considerably 
smaller at NLO and NNLO levels as compared with both the NLO higher twist terms obtained 
in [Q] and and the results of analysis obtained with no F-cuts imposed on the BCDMS data (see 
Fig. 2). 

4. Conclusions 

In the paper [jl]] the Jacobi polynomial expansion method developed in [Q, [10|] was used to per- 
form analysis of <2 2 -evolution of DIS structure function F2 by fitting all the existing to date reliable 
fixed-target experimental data that satisfy the cut x > 0.25. Based on the results of fitting, the QCD 
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Figure 2: Comparison of the HTC parameter h^x) obtained at LO, NLO and NNLO for hydrogen data 
when no Y cuts imposed on the BCDMS data. 



coupling constant value at the normalization point was evaluated. Starting with the reanalysis of 
BCDMS data by cutting off the points with large systematic errors it was shown [IJ, ||] that the val- 
ues of a s (Mz) rise sharply with the cuts on systematics imposed. The values of a s (M|) obtained in 
various fits are in agreement with each other. An outcome is that quite a similar result for a s (M|) 
was obtained §\§ in the analysis performed over BCDMS (with the cuts on systematics) and the rest 
of the data, thus permitting us to fit available data altogether. 

It turns out that for Q 2 > 3 GeV 2 the formulae of pure perturbative QCD (i.e. twist-two 
approximation accompanied by the target mass corrections) are enough to achieve good agreement 
with all the data analyzed. The reference result is then found to be 

a,(M|) = 0.1 167 ±0.0008 (stat) ±0.0018 (syst) ± 0.0007 (norm) 

= 0. 1 167 ± 0.0021 (total exp. error) . (4. 1) 



Upon adding twist-four corrections, QCD (i.e. first two coefficients of Wilson expansion) and 
the data are shown to be consistent with each other already at Q 2 = 1 GeV 2 , where the Wilson 
expansion begins to be applicable. This way we obtain for the coupling constant at Z mass peak: 

a s {Ml) = 0.1 167 ± 0.0007 (stat) ± 0.0020 (syst) ± 0.0005 (norm) 

= 0. 1 167 ± 0.0022 (total exp. error) . (4.2) 



Note that the above values ( |4.l| ) and (4^2) are to some extent stable [13] under the application 
of the "frozen" [14] and analytic [|^] modifications of the strong coupling constant, which as a rule 
lead to similar results (see [|l6|]). 

Note also that our results (4.1) and (4.2) for a $ (M|) are in good agreement with the world 
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average value for the coupling constant presented in the review [17], 



a s (Mf) =0.1 184 ± 0.0007. 



(4.3) 



Concerning the contributions of higher twist corrections in the present work the well-known 
x-shape of the twist-four corrections while going from intermediate to large values of the Bjorken 
variable x is well reproduced. 

A.K. is indebted to organizers for the possibility to present the talk which this paper is based 

on. 
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